Nonadiabatic production of spinor condensates with a QUIC trap 
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Motivated by the recent experimental observation of multi-component spinor condensates via a 
time- dependent quadrupole-Ioffe-configuration trap (QUIC trap), we provide a general framework 
for the investigation of nonadiabatic Landau-Zener dynamics of a hyperfine spin, e.g., from an 
atomic magnetic dipole moment coupled to a weak time-dependent magnetic (B-) field. The spin 
flipped population distribution, or the so-called Majorona formula is expressed in terms of system 
parameters and experimental observables; thus, provides much needed insight into the underlying 
mechanism for the production of spinor condensates due to nonadiabatic level crossings. 
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I. INTRODUCTION 

Magnetic traps play an important role in the study of 
atomic Bose-Einstein condensates (BEC) 1]. In a typical 
static magnetic trap, individual atomic spin couples to 
the spatial dependent magnetic (B-) field because of Zee- 
man effect. When an atom moves in a region where the 
direction of the B-field changes slowly and the strength 
of the B-field is sufficiently strong, according to Born- 
Oppenheimer approximation 0, U the atomic spin can 
follow the B-ficld adiabatically and remain in the same 
trapped eigen-state of the interaction Hamiltonian rela- 
tive to the instantaneous direction of the magnetic field 
B(r), where f is the center of mass position of the atom 
(or more precisely, that of the valence electron). In this 
case the atomic center of the mass experiences an effec- 
tive spatially-varying potential that is equal to the Zee- 
man energy and proportional to the strength of the B- 
field. 

For weak B-fields, when the atomic Zeeman energy is 
comparable to or less than the frequency of the direc- 
tional variation of the B-field felt by the moving atom, 
adiabatic dynamics cannot be followed anymore. As a 
result, nonadiabatic (Majorona) transitions 4] for the 
atomic spin may occur. Two potentially damaging ef- 
fects can cause nonadiabatic transitions. The first hap- 
pens when an atom enters a weak B-field region due 
to its translational motion in space. For instance, in a 
quadrupole trap, atoms in the weak field seeking state 
are accelerated towards the center of the trap where the 
B-ficld vanishes. Nonadiabatic transitions always occur 
in the vicinity of a zero B-field. To avoid this region of 
vanishing B-ficld or a spatial "hole," a number of meth- 
ods have been developed to effectively plug it, e.g., with 
the use of a far-off-resonant optical potential as an "op- 
tical plug" [| or the more famous time averaged orbiting 
potential (TOP) trap |fj. The second reason for nonadi- 
abatic transitions is due to the explicit time dependence 
of the B-field. Obviously nonadiabatic transitions may 
occur if the B-field changes rapidly with time. 

Recently, the atomic quantum gas group at Peking 



University (PKU) reported interesting observations of 
multi-component 87 Rb (F = 2) spinor condensates via 
switching off the B-fields of an initially spin polarized sin- 
gle component condensate in a QUIC trap 7] . The group 
of Prof. Chandra Raman at Georgia Tech also discovered 
counter-intuitive meta-stability when condensates were 
loaded into an "unplugged" magnetic quadruple trap . 
We decided to present our theoretical studies in the hope 
that the theoretical framework for spinor nonadiabatic 
level crossing dynamics may be of interest to other groups 
in the field of atomic quantum gases. In this paper, we 
will focus on the Peking University experiment in a time- 
dependent QUIC trap [!j ■ The more involved situation of 
a condensate in a quadruple trap will be discussed else- 
where |U| . According to the reported experiment [7J , the 
affected time-dependence for the B-field is relatively sim- 
ple. After a single component condensate was created in 
a QUIC trap, the various B-field generating currents were 
switched off in appropriately chosen orders. Whenever 
near-zero level-crossing occurs, multi-component spinor 
condensates are observed. 

This paper summarizes our treatment of level crossing 
dynamics for an atomic spin inside an external B-ficld. 
The theory is developed with respect to "the first sce- 
nario," where the vanishing B-field is due to the different 
time constants of decay for the B-fields from the QUIC 
coil and the bias coil after being shut off as discussed in 
Sec. II. An alternative scenario where the B-field zero is 
due to different time constants of the decaying B-fields 
from the quadruple coil and the Ioffe coil will be dis- 
cussed in Sec. III. Finally we conclude and provide a 
brief summary in Sec. IV. 

II. THE FIRST SCENARIO 

The magnetic trap used in their experiment Q is made 
up from two separate coils, a QUIC coil and a bias coil. 
The QUIC coil consists of a quadruple trap coil and an 
Ioffe coil in series as in the original QUIC trap 9J. The 
compensating coils for the earth's B-field are separate 
and always left on; thus will not be included explicitly 
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Quadrupole coils 




in our model. Before switching off, the magnetic B-field 
B^(r) created by the QUIC coil has the familiar config- 
uration of a Ioffe-Pritchard trap and can be expressed as 

B^(r) = B^(r)e ± + B^(r)e z , (1) 

where the axial and radial QUIC B-field components are 
B ( ^(r) and B^(r), respectively; 

B?(r) = gQ(0 )+igy , 

B?(r) = (2) 

with the unit vector ej^ defined as 
(-xe x + ye y ) 

e± = / 9 9 • W 

\l x + y 

B^ 1 and B^" denote the respective spatial derivatives for 
the B-ficlds. The right handed coordinate system as in 
Figure HI is chosen such that B$(0) > and B®' > 0. 

The bias field B A is in the z direction. It is created 
by the bias coils and approximately constant near the 
trap center 0. Before switching off, it can be expressed 
as B A (r) = -B A e z satisfying B®(0) > B A > 0. If 
the bias field is switched off first and the two compo- 
nents of the QUIC field are simultaneously switched off 
after a time interval of t[ nt , then at time t, the QUIC 
field becomes e~ t ^ TQ B'~ i {r) and the bias field becomes 

g-(t+t tat )/T A ^A^ ) Le>j both th g QTJ IC field and th g biag 

field are assumed to decrease exponentially with decay 
time constants tq and ta, and the quadruple field and 
the Ioffe field are assumed to decay with the same time 
constant. Assuming t — as the instant for shutting 
off the QUIC field, the total time-dependent B-field then 
becomes 

B(r,t) = e~ t / TQ B ( l(r)e ± 

+e-'/ r «5 2 Q (r> z - e-( t+t ^/ TA B A e z . (4) 



B±(r) is proportional to (x 2 +y 2 ) 1 / 2 near the trap center 
or the origin; therefore, we have B^(r), B A ^> £>5(j) 
and \B®(r) - B A \ > B®(r). Before switching off the 
QUIC field, the z component of B(r, t) takes a positive 
value B^(r) — B A much larger than the initial value of 
the transverse field B±(r). At t = 0, all atomic spins ini- 
tially are polarized, thus are the eigen-state \Mp = 2} of 
F z , i.e., the z component of the atomic hyperfine spin F. 
If the QUIC field and the bias field are switched off simul- 
taneously, i.e., tint = and tq ~ ta, the direction of the 
total B-field B(r, t) does not change with time although 
the strength of B(r, t) decreases after the switching off 
process. Nonadiabatic transitions do not occur in this 
case and the initial single component condensate remains 
a single component one. If the QUIC field and the bias 
field decrease with different time constants tq ta, the 

direction of B(r, t) changes with time and nonadiabatic 
level crossing arises. 

In the calculations to follow, we will make a simple 
approximation that the atomic spatial position does not 
change during the switching-off process. This allows for 
an easy calculation of nonadiabatic transition probabili- 
ties between different atomic spin states at a fixed spa- 
tial position r. This is well justified for the experiment 
of PKU, where level crossing occurs over a time window 
of ~ 10 2 (yus), during which an condensed atom moves a 
distance less than 0.1 (/im), provided its kinetic energy is 
-10 4 (Hz). 

As mentioned above, the z component of the B-field 
initially takes a large positive value. After the switching- 
off, the bias field decreases much slower than the QUIC 
field 0, i.e. we have ta 3> tq. At certain instant to, the 
value of e~ t / TQ i?P equals e ~^ t+tint ^ TA B A , which causes 
the z component of the total B-field to become zero. As 
a result of this, transitions from the state \Mp = 2) to 
other eigen-states of F z occur because of the finite trans- 
verse B-field e^ t / TQ B ( ^ in the vicinity of to. After to, 
the z component of the B-field becomes negative because 
e -(t+tint)/T A j e -t/T Q ^ 2 for a large enough t; the ab- 
solute value of the z component of the B-field can be- 
come again much larger than the transverse components 
of B(r,t) for t 3> to- Therefore, the probabilities for an 
atom in different eigen-states of F z can again take con- 
stant values in the long time limit. 

To compute the nonadiabatic level crossing rates, we 
note that transitions mainly occur in the near zero B-field 
region, i.e., for weak B-field. Thus, we only need to con- 
sider the linear Zeeman coupling of an atomic hyperfine 
spin. Our model Hamiltonian takes the simple form 

H = g Ff i B B{f, t) ■ F. (5) 

Here gp is the Lande g factor and [1b is the Bohr mag- 
neton. For 87 Rb atoms under consideration here, the 
spinor degree of freedom refers to the F = 2 manifold 
with gp = 1/2. In their experiment 0, the initial condi- 
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tion corresponds to 



|*(0)) = |M F = 2>. 



(6) 



At large t — * +00, the wave function can be expanded as 
2 

|*(t-+oo))= Yl C MF {r)e-^ M ^ t] \M F ), (7) 

AJ>=-2 

in the complete basis of F z along the initial quantiza- 
tion axis. Our problem is to find the steady population 
distribution PM F (r) = \Cm f {t)\ 2 in the long time limit. 

We will make use of the method of Hioe 11] to calcu- 
late the finial state probability distribution due to nona- 
diabatic level crossing of a high spin. Because of the 
rotational symmetry of our model system JSJ), it can be 
mapped onto a spin 1/2 spinor with the same type of 
coupling, described by a Hamiltonian 



(8) 



where a is the familiar spin 1/2 Pauli matrix vector. The 
initial condition for the spin 1/2 state is 



lv(0)> = [i,o] T , 

and the finial state can be denoted as 



\ip{t — > +00)) = [a(r)e 



/3(r)e 



(9) 



(10) 



Upon solving this two state problem, Pm f (r) can be 
found easily according to the rotation group represen- 
tation elements as in Hioe jn|. Apart from a globe 
phase factor, the evolution operator corresponding to the 
Hamiltonian JSJ) can be expressed as — exp[—in-F9]; 
while the one corresponding to the Hamiltonian (JHJ) is 
jji 1 / 2 ) = exp[— in ■ (a/ 2) 9]. The unit vector n and the 
angle 9 are determined by B(f, t). Therefore, Z^ 1 / 2 ) and 
Z)( 2 ) are the representation matrixes (D matrixes) of the 
same rotation operation. The transition probabilities 



\D 



(2) 



|2 _ 1 r^ 1 / 2 ) 1 2 According to the representa- 



,(2) 



and \D 



(2) 



Pm f and |a(r)| can be rewritten as Pm f 
and \a{r)\ 2 = \D^ /21/2l 

tion theory of SO(3) group 12], \D^ 2 1 ^ 2 , \^ Mf 2 
are functions of sin[/3/2] and cos[/3/2] with (3 one of the 
three Euler angles of the rotation. Although we do not 
know the values of n, 9, and (3, we can express the tran- 
sition probability Pm f in terms of A(r) = |a(r)| 2 ; as, 

P 2 (r) = A(f)\ 

P 1 (r) = 4A(rf[l - A(r)}, 

P (r) - 6A(r) 2 [l-A(r)} 2 , 

P-i(r) = AA(r)[l-A(f)} 3 , 

P_ 2 (r) = [l-A(r)]\ 



(11) 



The two state problem can be solved accurately with 
the Landau-Zener formula |l3[ . To this end, we reexpress 
the Hamiltonian JBl as 



H^t] = q[t\K[t], 



(12) 



with the "normalized" Hamiltonian 

h a [t] = gj_crj_ + (g 
and the parameters 



(13) 



9± = ^gF^B B ±:{r), 

g? - IgpmBfir), 

£ = (ta - tq)/t a t q . 
We define a new time variable 

s = -T Q e-* /r «; 
then, the time-dependent Schrodinger equation 

id t \ l p[t])=H lT {t]\ip.[t}), 

becomes 

id s \ip{s)) = h(s)\ip(s)), 



(14) 
(15) 
(16) 
(17) 



with <p(s)) = |(^[t(s)]) and h(s) = h[t(s)]. Consequently, 
the time interval of the dynamics t 6 [0, 00] is mapped 
into s E [— tq, 0]. 

As stated above, the z component of B(r, t) takes large 
positive and negative values, respectively, at t = and 
t = 00. Therefore, at s = — tq and s = 0, the condition 



> <7±, 



(18) 



is satisfied while g® — e~^ s ^g^ takes positive and neg- 
ative values, respectively. In the Landau-Zener approx- 
imation, a linear approximation is always assumed for 
the different energy levels. We find the value sq at the 
crossing point t is given by 



«o = -TQq(t ) = -tq ( % 

\g? 



At s = so when the longitudinal B-ficld vanishes 



t 91 



-£*(«o)„A 



0, 



(19) 



(20) 



a linear approximation to the energy levels simply leads 
to 



(21) 



with 



<g? 



Using the Landau-Zener formula, we immediately find 



A{f) = exp — 7r 



I9x| 



(22) 
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For ta > tq, we find £tq w 1, which then leads to 



A(f, t h 



exp 



= exp 



\g±\ 2 gtr Q 



-TTg F fJ, B B^(f)TQ 



Bj 2 (f) _ 



2B? 2 (r) 



— . (23) B 



_10 

t 

o 

i 8 



Obviously for a large enough time interval tint such that 
e -ttat/TA < i 5 we have A(r, t in t) « 1 and P 2 (r) « 1. A 
single component condensate remains a single component 
one. In fact, if e~* int / TA <C 1, the bias field has already 
decreased to zero when the QUIC field is switched off. 
Thus, during the switching-off of the QUIC field, the 
direction of the B-field does not change and nonadiabatic 
transitions cannot occur. 

In the experiment of PKU, 7], the various trap pa- 
rameters take the following values: B®(0) = 9 (Gauss), 
#A = 7 45 (Gauss), Bf = 4.9 x 10 2 (Gauss-cm- 2 ), 

= 3.0 x 10 2 (Gauss-cm" 1 ), r Q = 40 (/is), r A = 3 
(ms). Before switching-off, the center of the QUIC trap 
is at r*o = (0,5/im, 0). Substituting the above coefficient 
A(f) of Eq. (|23|l into Eq. Qllfl. we arrive at a simple esti- 
mate for the population distribution Pm f (^b)- More pre- 
cisely, we can estimate the population distribution with 
the following, 



N 



Mr 



Pm f (r)p(r)dr, 



(24) 



where p(r) is the density profile of the trapped gas cloud. 
Figure [3 shows the typical dependence of such result on 
the time interval t- ln f In the calculation of Figure 
we set p(r) to be the atomic density distribution given 
by the Thomas-Fermi approximation corresponding to 
the initial values of the QUIC field and the bias field. 
Namely, the spatial motion of the atoms is omitted. This 
approximation is based on the fact that the decay time 
(3ms) is shorter than the period (~4. 5-7. 5ms) of the trap 
potential and seems to be a bit crude. To obtain a more 
accurate estimation of the atomic population, variations 
of the atomic spatial distribution in the decay process of 
the bias field should be considered fully . 

The above result is based on the approximation that 
atoms do not move during the switching-off process. If 
atomic motion is included, more accurate population 
distribution can be calculated by solving the multiple- 
component Gross-Pitaevskii equation including the time- 
dependent B-field. A detailed comparison of these two 
approaches is given in Ref. [Tojj. Overall, we find the 
approximate Landau-Zener solution discussed here holds 
well for the parameter regimes of the experiment Q . 



III. AN ALTERNATIVE SCENARIO 

The QUIC coil consists of a pair of quadruple coils 
and an Ioffe coil. The B-field B® of the QUIC trap is 
the sum of the B-fields B qd from the quadruple coils and 
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FIG. 2. (Color online) A typical dependence 

of the population distributions on the time 
ti nt . The unit of the atomic number is 10~ 4 . 

B u from the Ioffe coil. In the previous section, we simply 
assumed the B-fields generated by these two sets of coils 
decrease synchronously after shutting off electric currents 
. However, as was discovered in the experiment 0, the 
magnetic fields B qd and B li do not always decay with 
the same time constant despite the fact that the two sets 
of coils forming the QUIC trap are in series. Assuming 
different time constants for the decay of B qd and B , the 
finial population distribution needs to be re-calculated. 

Before the QUIC field is switched off, the components 
of the B-fields B qd and B li are functions of the atomic 
position, explicitly given by 

Bf(r) 
B*(r) 
Bf(r) 



B q 



(1/ 



B*<f) 
Bf{r) 



zo), 

-2B qd 'x, 

2B qd ')a 



z )+B?(0), 



(25) 



B qd 'y, 

where zq is the distance between the center of the QUIC 
trap (in the absence of gravity) and the center of the 
quadruple trap. 

In their experiment Q, B qdl is about 150 (Gauss- 
cm -1 ) and zq is 0.75 (cm). Therefore, in the region 
near the center of the QUIC trap, we have B qd = —107 
(Gauss). From B$(Q) = 9 (Gauss), we find B l J = 116 
(Gauss). 

If B u and B qd decrease with different time constants 
th and r q d after switching-off, the B-field from the 
quadruple coils becomes 

(Bfe z + Bfe y + Bf e^e"*/^ (26) 

and the B-field created by the Ioffe coil becomes 

(B^ + B^ey + B^e-^. 



(27) 
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In this section, we assume the time interval t; n t be- 
tween the switching-off of the B-fields B® and B A is suf- 
ficiently long, i.e., when the QUIC field is switched off, 
the bias field has already decreased to zero. Since the 
B-fields Bf and B l J have different signs, at a time to 
the condition 



£lfg-io/Trf _|_ ^qd e -io/r q d _ q 



(28) 



can be satisfied and the z component of the total B- 
field becomes zero. As before, nonadiabatic transitions 
happen mainly in the temporal region near t . Assuming 



Tlf 



Tqd 

2 : 



(29) 



the nonadiabatic transition probability A in the spin- 
1/2 case can again be calculated with the Landau-Zcncr 
method used previously provided that to < t^. In the 
present case, we find to < 0.6rif. Thus, we obtain 



A(r) = exp 



9 Ff iBruEi«\Bl f \ + Bj ! \Bf\) 
IOT 



,(30) 



where I = x,y, the counterpart of the parameter A in the 
first scenario. 

Now we discuss a special case. We assume the bias field 
is switched off adiabatically, such that the atomic cloud 
follows the variation of the total B-field, and moves to the 
region near the center fi of the QUIC trap (in the ab- 
sence of the bias field). Since the Landau-Zener method 
is based on the assumption that the atoms are located in 
the region where the QUIC field lies approximately along 
the z axis before being switched off, the factor A given 
in Eq. Q30f) is applicable if the trap center r\ is near the 
z axis so that |_B^,^j| is much smaller than \B®\. For 

practical values of B and J5 qd , the above condition is 
satisfied and a good estimate for the transition probabil- 
ity can again be given by Eq. (|30|l . 

From the directions of the electric currents in the 
quadruple coils and the Ioffe coil as shown in Figure Q 
the B-fields By 1 and B^j d are found to have the same sign 
while the fields f?* f and B^ d have opposite signs. There- 
fore, after switching off, B^ decays much slower than 

B®. At time to when B® = 0, B® has the same order of 
magnitude as its initial value. Therefore, if in the region 
near r*i, B® is sufficiently large, the direction of B® may 
be changed very slowly during the switching-off process 
of B® so that the atomic spin state can be adiabatically 



flipped. Then, we have A rj and P-2 ~ 1. For in- 
stance, in the experiment of Ref. 0, = (0, 30^m, 0) 
and ~ 0.45(Gauss). In this case 85 percent of the 
atoms can be switched to the state \Mp = —2). 

A second case of some interest is when the bias field 
is switched off suddenly. Once the bias field is turned 
off, the atoms begin to oscillate in the new QUIC trap 
centered at r\ . If the time interval between the switching 
off of B A and B® is tint , then the population distribution 
can be estimated as 



N Mf = / drp{r,t int )P MF [Mr)}- 



(31) 



Here, p(f, t int ) is the density distribution of atoms in the 
QUIC trap at the time when the QUIC field is switched 
off. 



IV. CONCLUSIONS 

In conclusion, we have presented a detailed theoretical 
treatment for the nonadiabatic level crossing dynamics 
of an atomic spin coupled to a time dependent magnetic 
field. When applied to the condensate experiments in 
a modified QUIC trap 7], our theory provides a sat- 
isfactory explanation for the observed multi-component 
spinor condensates when the trapping B-fields were shut- 
off. In the broad context of condensate wave function 
engineering and atom optics with degenerate quantum 
gases, our work provides useful insights for experiments. 
For example, in some proposals [14| and experiments 
on the creation of vortex states in a condensate, the in- 
ternal atomic hyperfine state is slated to adiabatically 
follow the external magnetic field and be changed from 
\^f)z to | — mp) z . Our method can then also be used to 
estimate the nonadiabatic effects in these proposals and 
experiments [T3. IT5l |. 
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